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Abstract 

By using a relativistic framework and accurate nuclear spectral func- 
tions we evaluate the ratio F^A/AF^iy of deep inelastic neutrino scattering. 
Parametrizations of this ratio for different values of are provided. These 
results should be useful for taking into account the nuclear effects in analyses 
of experimental data in neutrino reactions in nuclear targets, and test QCD 
predictions for the nucleon structure functions. In particular, the size of the 
nuclear corrections is of the same order of magnitude as the size of the QCD 
corrections to the Gross-Llewellyn Smith sum rule. 
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Typeset using REVT^ 



I. INTRODUCTION 



The deep inelastic scattering (DIS) experiments done with neutrino beams can now 
provide data with an accuracy comparable to the analogous processes studied with charged 
leptons [HJ^. Due to the low cross sections, for neutrino interaction with matter, these data 
are collected mostly from experiments done on heavy nuclei like iron 0. 

Since the discovery of the EMC effect the deviation of the nuclear structure function 
F2{x) with respect to the deuteron, many theoretical attempts were made to explain it (for 



a review see In the conventional nuclear physics models, the usual effects considered 

are Fermi motion 0] and binding of the nucleons in the nuclei and the contribution of 
the virtual meson cloud [P,p!0[]. 

For the neutrino processes, the study of possible nuclear effects has received much less 
attention, probably because of the lack of accurate measurements in the past. In present 
analyses of neutrino experiments, the high x data [x > 0.7), where Fermi motion effects are 
expected to be important ||^, are generally avoided, especially in the case of nuclear data 
0]. However, it has been shown that even in the case of the deuteron, the Fermi motion 
effects could be important [|1T|. It is interesting to note that the older fits of Diemoz et al. 
|, show that even for lower x (x ~ 0.5), the nuclear effects could be important. Most of 



the analyses done for the neutrino data assume the same nuclear corrections as applicable 
to the charged lepton experiments in their global fits of the parton distributions |jl3|,|14[, 
while allowing for an overall normalisation of the data. This is not appropriate for the 
parity violating structure function F^j^{x), in the neutrino scattering data, which has no 
counterpart in the charged lepton case. Even in the case of the structure function Fr 



2A 



X 



which is related to the F2a{x), of the EMC effect, the nuclear corrections would not be the 
same due to the different contribution of the sea quarks [|15 



Recently, there have been some calculations of the nuclear effects in the F^^{x) structure 
function [lT6| , p!7| . It is found by Sidorov and Tokarev (ST) [|l6l that the nuclear effects in 
are quite small in the case of the deuteron, a result similar to the one found in charged 



[X 



lepton scattering in -F2d(^) [HI) 19 1- Keeping in view the result that the nuclear effects in 



F2{x) are quite different for the deuteron and heavy nuclei, the assumption that they are 
the same for Fn( done by ST is therefore not justified. 

In the calculation of Kulagin |[T^, the iron structure function is evaluated and sizeable 
deviations from the deuteron one were found. The effect on the Gross-Llewellyn Smith 
(GLS) sum rule [20| is on the contrary very small. This is so because in Kulagin's approach, 
the spinors are normalized in order not to affect this sum rule. It is similar to the prescription 
introduced by means of the so called flux factor in order to take into account relativistic 
effects in non relativistic formalisms. It is interesting to note that in the fit that the CCFR 
collaboration made to the xF2{x) and xFslx) neutrino data P,p2|, the integral of the valence 
quark distributions is not 3, p2]. This apparent lack of normalization of the ^3^1(0;) in the 
nuclear data is of crucial importance in order to extract the strong coupling constant, a^, 
from the analysis of the GLS sum rule, as done in P3 . 



Recently, a new calculation of the nuclear effects in the case of the charged lepton scatter- 
ing was made |ll0| , |2^ . In this calculation a relativistic many-body approach for the nucleon 



contribution was developed in order to avoid the use of the flux factor mentioned before ^1 
The meson cloud contribution was found to give relatively important contributions. Good 
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agreement was found with the experimental data. For the structure function Fsa{x), the 
virtual meson cloud gives no contribution, since it is related to the difference of the quark 
and antiquark distributions. In the present paper we calculate the nuclear effects to ^3^(0:) 
in the same formalism of ||10|,^ . 

The structure of the paper is as follows. We briefly present the formalism for deep 
inelastic neutrino scattering in section II. In section III we show how nuclear effects are 
calculated. In section IV we present the results and state our conclusions. 



II. DEEP INELASTIC NEUTRINO-NUCLEON SCATTERING 

The invariant matrix element for the Feynman diagram shown in Fig. 1 corresponding 
to the inelastic neutrino-nucleon scattering is written as 

-zT= (^^^ - l.)uAk) (^i^l^) (^l-^-l^) • (1) 

The cross section for the neutrino scattering, o"'^, is then given by 

1 2my 2M r d^k' 2m 

a 



Vrei 2E,{k) 2E{p) J (27r)3 2E^{k') 



]eb \ J, 



\J (2vr)3|^i \2E[ 

N 

(2vr)V(p + fc-fc'-Ep:), (2) 



i=l 

where / stands for fermions and h for bosons in the final state X. The index i is split in Z, j 
for fermions and bosons respectively. 

By inserting expression (|I]) in Eq. (^), summing over the final spins and averaging over 
the initial spins, the cross section can be expressed in the rest frame of the nucleon, with 

E' referring to the outgoing muon, as 

dQ'dE'-{2ny \k\ W-ml.) + ^5 J , (3) 

where the lepton tensors, L"^ and Lg^ are given by 

L^P = k'^k''^ + k^k"" - {kk')g"'^ , 

Lf = -te'^^P-kX , (4) 
and the hadronic tensor is defined as 

N 

X {X\J^\N){X\MNn2nr6\p + g - > (5) 

i=l 
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where q is the momentum of the virtual W , Sp the spin of the nucleon and Si the spin of the 
fermions in X. 

In the case of antineutrino scattering, the invariant T matrix is given by 
and the cross section reads 



dn'dE' (27r)2 |fc| \q'^-m 



w . 



The usual convention is to express the hadronic tensor as ||25 



+ {Paqfs + qaPf3)W^^''^ + — {p^qp - qaVpW^'^''^ , (8) 

where W^'^'^^ are the structure functions, which depend on the scalars and pq. The terms 
depending on W4, and Wq in Eq. (||) do not contribute to the cross section in the DIS 



Defining the variables 



Q'' = -q^] 1^=^; a; = -^; y = , (9) 



we can write in the Bjorken limit 

^ - ^ ^ (- - 1:^) -^'"'<^' -^<- ^/^)-^^'''<^ ' 

(10) 

where the + (— ) sign stands for the neutrino (antineutrino) cross section, and the F^^^\x) 
are dimensionless structure functions defined as 

Ff(') = Miyf(^); F!^^''^ = vW!^^'^- F^^'^ = vW^^'\ (11) 

In the quark parton model, the structure functions are expressed in terms of the quark 
distributions, q{x) and q{x). Using Callan-Gross relation, these structure functions are 
given as 

2xF'(^ = 2xF,\ = F^^ = F^j, = x{qix) + qix)) , 
= Qi^) ~ Qi^) + 2s (x) - 2c(x) , 
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F!j, = q{x)-q{x)-2s{x) + 2c{x), 
Fs{x) = ^[F^n{x) + F^j,{x)] = q{x) - q{x) = u,{x) + d,{x) . (12) 



Here 



q{x) = u{x) + d{x) + s{x) + c{x) 

q{x) = u{x) + d{x) + s{x) + c{x) . (13) 

The neutrino structure function Fl'jy is related with the analogous one for charged lepton 
scattering by 

^ ^ ^ / _ 3 s{x) + s{x) - c{x) - c{x) \ 

F^M 18 1^ 5 q{x) + q{x) J ' ^ ^ 

In the deep inelastic region, quark distributions satisfy sum rules in order to give the correct 
charge, strangeness and charm of the proton and the neutron. This implies the constraint 

1 

F3{x)dx = 3, (15) 







known as the Gross-Llewellyn Smith (GSL) sum rule EC 



QCD corrections modify the expression of the structure functions in terms of the quark 
distributions, Eqs. (P^). At leading order, Eqs. (|T2|) retain the same structure but the quark 
distributions evolve with according to the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi) evolution equations When higher orders are considered, Eqs. (12) are modified 
including explicitly the strong coupling constant, asiQ"^)- 

The GLS sum rule is modified and in the NLO approximation is given by 

Fs{x) dx = 3{1 - ^^^) . (16) 
Jo vr 



Further higher order corrections up to order a|(Q ) can be found in |27] 



III. NUCLEAR EFFECTS IN NEUTRINO SCATTERING 

In order to calculate the neutrino-nucleus cross section we first evaluate the related 
neutrino self-energy in the medium. Fig. 2. The self-energy is given by 

G f d^q _ , . ^' + m^ 



H) ^"-'{q) , (17) 



g2 - 



where 11"^ (g) is the W self-energy in the medium. 
Eq. ([T7| ) can be rewritten in the form 
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The probability per unit time for the neutrino to colhde with nucleons when traveUing 



through nuclear matter is |10 



2777 

r(fc) = --^ImS(fc), (19) 

and the cross section for an element of volume d^r in the nucleus is 

da = VdtdS = V— dldS = — d^r = 
dl V 

^ r d^r = Im S d^r. (20) 

|fc| \k\ 

Im S can be easily evaluated from Eq. ([18|) by means of the Cutkosky rules |^ 



S(A;) ^ 2i Im S(A;) 

D{k') — > 2i6'(A;'°) Im _D(fc') (boson propagator) , . 

W^iq) 2i9lq°)lmW{q) ^ ^ 

G{p) 2iO{jP) Im G{p) (fermion propagator) 



and we get 

rfV^ G 4 / mw 



dVL'dE' V2(27r)3 \^ 



2 
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W 



d\ Im n,^(g) . (22) 



Comparing Eq. (||) and Eq. (|22|) we see that 

V2 _ 

Next, we evaluate the W self-energy in the medium 



W^/(?) = -— TtV Idhluvir^iq). (23) 
TT Grriw J 



--^'"-(-'/|f,.«(P)i:En/|| 



niG,(p;)n iDAp'i) (-^) (A'|J'1A'>(.Y|J''|/V>- 



X 

' I j 

\4r4/„ I „ I 



x(27r)V(g + p-Sf=iK), (24) 

where we have a minus sign because of the necessary fermion loop. In the antineutrino case 
the expressions obtained are very similar. V^^ appears, as in Eq. (|^), with a minus sign in 
front and in the W self-energy, Eq. (H), we have (X|J^|A^), instead of {X\Ja\N). From 
now on we will always speak of the average of neutrino and antineutrino structure functions 
and will omit the superscripts v and v. For the nucleon propagator in the medium, G{p), 
we take a relativistic version |]ID|, which can be written as 
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M 



+ 



a uo— 



UJ — 17] 



duo- 



p'^ — UJ + ii] 



(25) 



Sh{uj,p) and Sp{uj,p) being the hole and particle spectral functions respectively, which are 
taken from the work of [|10],^|. 



We use the local density approximation in which the spectral functions depend on the 
density of the point of the nucleus at which they are evaluated. In our formalism we use 
spectral functions for symmetric nuclear matter. The normalization of the hole spectral 
function is given by 



d^p 
(2^ 



S'/i(u;,p, kpir)) du = A. 



(26) 



where kp{r) = [3'7r^p(r)/2]^/'^ is the local Fermi momentum at the point r. All our calcula- 
tions are done for ^^Fe. The density for this nucleus is expressed as a two Fermi parameter 
distribution given in |^ . 

We now calculate Im 11"^ (g) by using Cutkosky rules, Eq. (|2TD , the expression of Eq. (^5|) 
for the nucleon propagator in the medium, free propagators for particles in the final state 
and by means of Eq. (^31) we have the hadronic tensor 



Wf 



4 / d^r 



d^p M 



(27r)3 E{p) 



dp'SH{p',p)Wf{p,q). 



(27) 



In order to evaluate F^a, we calculate the components xy on both sides of Eq. (p7|). We 
have, using €0123 = 1, and taking q along the z axis, as usual. 



W7 



2M 



-QzWsA 



A 



and for the right hand side we will have for the moving nucleon 



^xy _ PxPy I 



Since we have 



M2 ^^^^^'-^^ 2M2 

W3N{p,q) = F3n{xn) 



(28) 



(29) 



pq 
M 



(30) 



with X as defined in Eq. (^ and xat is the Bjorken variable expressed in terms of the nucleon 
variables , (p°,p), in the nucleus 



Xn 



2pq 



(31) 



we obtain the expression for F^^ix) in the Bjorken limit 
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A 



4 / d-^r 



d^p M 



(27r)3 E{p) 



' dp'S,{p',p)^ 

-oo X 



PoQz - PzQo 



(32) 



where the contribution of W2 appearing in Eq. (|2^) vanishes after momentum integration. 
Defining 7 as 



(33) 



we get 



A 



4 / d^r 



d^p M 



(27r)3 E{p) J-00 



dp'Sr,{p\p) 



Pol - Pz 

,{po-Pzh, 



FMxn) (34) 



IV. RESULTS AND CONCLUSION 



We begin by showing the results for = F^^/AF^i^f and comparing them with the corre- 
sponding ratio R2 = F^j^/AF^j^ in Fig. 3. For F^^ and F^j^ we have taken the parametriza- 



tion of ||31fl. For the meson structure functions that contribute to R2 we have used the 



parametrization given in ^2,\. We can see some similarities in the region around x = 0.5 ~ 0.6 



where the ratio is smaller than unity and which is mostly due to nuclear binding effects, 



as discussed in detail in [p|,pl|. Similarities in the region of x > 0.6, where the ratio shows 
a fast increase, are also apparent, and they are mostly due to the effect of Fermi motion 
||7|JT0[|. The differences between the neutrino and charged lepton ratios are more important 



at values of x < 0.6. These differences are mostly due to the lack of meson renormalization 
effects in the neutrino structure function, as can be seen in the figure, where we also show 
the results for R2 without the meson renormalization effects. 

We should note that the nuclear effects are sizeable, with values of R3 around 0.8 in 
the region of x ~ 0.6 and around 0.9 for low values of x. These nuclear corrections are 



considerably larger than those found for the deuteron [0, as it was also the case for R2 in 
the charged lepton case. 

In Fig. 4 we show results for R3 for different values of Q^. We observe that for low 
values of x the results for R3 are rather independent of Q^, but this is not the case at large 
values of x, where there are substantial differences. Given the fact that the most important 
contribution to the GLS sum rule comes from the values of x < 0.4, the results of the 
figure indicate that nuclear effects can reduce the GLS sum by about 10%. Since one of 
the purposes of this work is to facilitate the task of experimentalists in analyses of the GLS 
sum rule and other QCD predictions for -F37V, we provide here an easy parametrization of 
i?3 which can serve to induce the results for F^jy from the measured nuclear data of ^3^^. In 
table I we give the parameters of a fit of the ratio -R3(x) for the values of 5, 30 and 50 
GeV^. The form used is R^ix) = A{B — x)°/(l — x)^. For larger values of results are 
similar to those for 50 GeV^. 

We also want to stress that the factor {poj — Pz)/ {Po —Pz)l appearing in Eq. ( P^D which 
goes to unity in the Bjorken limit (see Eq. (0)), is not negligible for values of ~ 4M^ 
and produces changes of around 10% in R^ for values of x ~ 0.8 — 1 at = 5 GeV^. 



8 



In Fig. 5 we compare our results with those of Kulagin and ST [0. It is quite clear 
that both for the case of Kulagin and the present results the differences with the deuteron 
results of ST are significant. Our results and those of Kulagin are qualitatively similar 
around x ~ 0.4 — 0.5, but they divert from each other both at large x and small x. The 



discrepancies at large x are not surprising since it was shown in that the results for 
F2A at large x were very sensitive to the parametrization used for F2N, which is different in 
and the present case. These differences, however, should not be very relevant for tests 
of the GLS sum rule, since the contribution to the sum rule from the region of a; > 0.6 is 
very small. The discrepancies for values of x < 0.4 are more relevant in connection with 
this sum rule since 10% effects, typical differences at low x between both approaches, are 
of the same order of magnitude as the QCD corrections to the GLS. The main reason for 
the discrepancies in the two results should be seen from the different approaches followed. 
In ^ the spinors are normalized in a way to force baryon number conservation in a non 
relativistic formalism, in line of the idea of the flux factor of |^ to introduce relativistic 
corrections in a nonrelativistic formalism. The advantage of the present work is that it 
uses a relativistic formalism from the beginning, in which baryon number conservation is 
automatically fulfilled. The differences of 10% in the region of small x between the two 
aproaches is about the same one as found between our approach and those using the flux 



factor in studies of |[T0|JT8 



In summary, we present here results for based on a relativistic approach and the 
use of an accurate nuclear spectral function which reproduces nuclear bindings, momentum 
distributions, etc. Within this formalism we were able to reproduce the EMC effect for 
different nuclei and the experimental results of F2A in the region of x > 1. We have used 
here the same formalism to evaluate ^3^. There are some differences with respect to F2A 
in the nuclear effects, mostly due to the absence of meson renormalization corrections in 
F^A- The results obtained are parametrized in an easy way to facilitate future analyses of 
experimentalists in order to deduce the elementary F^n structure function from the nuclear 
measurements and be able to test QCD corrections and obtain reliable values of as- 
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TABLES 



TABLE 1. Fits of the ratio R^ix.Q'^) = F3a{x,Q'^)/AF3n{x,Q'^) for different values of 



Tli(^ form of Iho [il is Ry,(,r) = A(D - x)" /{I - xY 



Q2 (GeV2) 


A 


B 


a 


13 


5 


0.691 


1.29 


1.21 


0.688 


30 


0.647 


1.06 


6.88 


6.19 


50 


0.599 


1.07 


6.43 


5.64 
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FIGURES 



FIG. 1. Diagram for the process of deep inelastic neutrino- nucleon scattering. 

FIG. 2. Self-energy diagram of the neutrino in the nuclear medium associated with the process 
of deep inelastic neutrino-nucleon scattering. 

FIG. 3. Solid line: results for the ratio F3A{x)/AF3j\f(x); dashed line: results for the ratio 
F2j^{x) / AF2j\f{x) including the contribution of the nucleons and the mesons; dotted-dashed line: 
results for the ratio Fj^(a;)/AF2Ar(^) including only the contribution of the nucleons. All curves 
are evaluated at = 50 GeV^. 

FIG. 4. Results for the ratio Fsa{x) / AF^n {x) at different values of Q^. Dotted-dashed line: 
ratio at = 50 GeV^; long dashed hue: ratio at = 30 GeV^; sohd line: ratio at = 5 GeV^; 



short dashed line: ratio at Q = 5 GeV but setting 7 equal to 1 in Eq.(34) 



FIG. 5. Results for the ratio F^a^x) / AF^^ {x) at = 5 GeV^ by different authors. Sohd line: 
this work; dashed line: Kulagin [|l^; dotted-dashed line: Sidorov and Tokarev |0] (for deuteron). 
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(x) H 



